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Abstract 



| We consider the monomer-dimer model on sequences of random 

graphs locally convergent to trees. We prove that the monomer density 
converges almost surely, in the thermodynamic limit, to an analytic 
function of the monomer activity. We characterise this limit as the 
expectation of the solution of a fixed point distributional equation. 



1 Introduction 



Each way to fully cover the vertices of a finite graph G by non-overlapping 
monomers (molecules which occupy a single vertex) and dimers (molecules 
which occupy two adjacent vertices) is called a monomer-dimer configura- 
tion. Associating to each of those configurations a probability proportional 
\ to the product of a factor w' > for each dimer and a factor x' > for 

■ each monomer defines a monomer-dimer model on the graph. It is easily 

. seen that the monomer-dimer probability measure depends only on the fac- 

tor x = x'/vV. What one is mainly interested in are the monomer (and 
dimer) densities, i.e. the average number of monomers (dimers) per site. 
Monomer-dimer models were introduced in the last century in the physics 
^ literature to study the statistical mechanics problem of diatomic oxygen 

adsorption on tungsten pQ and similar phenomena (see [2] and references 
therein). Important rigorous results were obtained by Heilmann and Lieb in 
[31 E]) where in particular the absence of phase transition for the pressure as a 
function of x (and of the monomer density too) was proved for all positive x. 
Furthermore in |2j exact solutions were given for specific topologies like the 
one-dimensional (with free and periodic boundary conditions), the complete 
graph and the Bethe lattice. Previously exact solutions on two-dimensional 
lattices where found by Kasteleyn, Fisher and Temperley [H |6] for the 
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pure dimer problem, i.e. the problem of counting configurations with no 
monomers. 

In this paper we study the statistical mechanics properties of monomer- 
dimer systems on locally tree-like random graphs and compute their monomer 
density in the thermodynamic limit for all positive x. The class of diluted 
graphs that we cover is the same for which the exact solution of the ferro- 
magnetic Ising model was recently found by Dembo and Montanari [7\ using 
the local weak convergence strategy developed in [8]; precisely we consider 
random graphs (G n ) n& ^ locally convergent to a unimodular Galton- Watson 
tree T(P, p) and with finite second moment of the asymptotic degree distri- 
bution P. A remarkable example is the Erdos-Renyi graph, i.e. the complete 
graph randomly diluted with i.i.d. Bernoulli edges and average degree c. 

In the Erdos-Renyi case our main result is the proof that the monomer 
density e(x) exists in the thermodynamic limit and it turns out to be the 
expected value of a random variable Y x whose distribution is defined as the 
only fixed point supported in [0, 1] of the distributional equation 

Y 1 X - (1) 

where the (li)igN are i-i-d- copies of Y, K is Poisson(c)-distributed and 
independent of (Y^jgjj. A side-result, yet a crucial one, of our analysis is 
that the solution Y x is reached monotonically in the number of iterations 
of the previous equation. More precisely, starting from 1^ = 1, the even 
iterations decrease monotonically, the odd ones increase monotonically (see 
Fig. 2), their difference shrinks to zero and their common limit is an analytic 
function of x. 

Our results are built on the Heilmann-Lieb recursion relation for the 
partition function Zq{x) of a monomer-dimer system |2j. Given a finite 
graph G, a root vertex o and its neighbours i, it holds: 

Z G (x) = x Z G - {x) + ^2 ZG-o-i{x) ■ (2) 

In [9] it is shown how to rewrite the identity ([2]) in terms of the probability 
lZ x (G,o) of having a monomer in o. We use this form to deduce the distri- 
butional identity ([1]) for Y x := lim^oo 1Z x (T(r), o), where T(r) is a random 
tree with root o, r generations and i.i.d. Poisson(c) offspring sizes. Our re- 
sults rely on a correlation inequality method, which permits a local study 
of the monomer-dimer system and is also at the origin of the monotonic- 
ity property described before. Analytic continuation techniques are used in 
order to extend results from "large" x to all positive x. 

Our results extend those of Bordenave, Lelarge, Salez [S], which are valid 
for graphs with bounded degree and are generalised to arbitrary degree only 
for x — )■ 0, i.e. for the maximum matching problem (see also [TO] for this 



2 



problem in the Erdos-Renyi case). A complete theoretical physics picture of 
the monomer-dimer model on sparse random graphs was given by Zdeborova 
and Mezard in where several quantities were computed, including the 
pressure of the model, using the so called replica-symmetric version of the 
cavity method. Our work provides a rigorous proof of their picture, in 
particular it may be shown [12] that the cavity method pressure coincides 
with the one that can be obtained with our methods. 

The paper is organised as follows. Section 2 introduces the definitions 
and the basic properties of the monomer-dimer models, including their well- 
posedness with stability bounds for the pressure, the main recursion relation 
for lZ x (G,o), the analyticity property of its solutions, and the correlation 
inequalities for locally tree-like graphs at even and odd tree depth. Section 
3 studies the model on trees, in particular the solution on a Galton- Watson 
tree is found in Theorem [1] and its corollaries. Section 4 presents the general 
solution on locally tree- like graphs in Theorem [2] and its corollaries. Section 
5 displays lower and upper bounds for the monomer density in the Erdos- 
Renyi case, obtained by iterating the recursion relation ([1]) an odd and even 
number of times. 

2 Definitions and general properties of the monomer- 
dimer model 

Let G = (V, E) be a finite simple graph with vertex set V and edge set 
E C {uv = {u, v}\u,v £ V, u ^ v} . 

Definition 1. A dimeric configuration on the graph G is a family of edges 
D C E no two of which have a vertex in common. Given a dimeric configu- 
ration D, the associated monomeric configuration is the set of free vertices: 

jg G {D) := {v G V\ VuGVuv i D}. 

We say that the edges in the dimeric configuration D are occupied by a 
dimer, while the vertices in the monomeric configuration ^Kq(D) are occu- 
pied by a monomer. Notice that \^q(D)\ = \ V\ — 2 \D\ . 



Fig. 1: The bold edges in the left figure form a dimeric configuration on the graph, 
while those in the right figure do not. 
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Definition 2. Let $>g be the set of all possible dimeric configurations on 
the graph G . For a given value of the parameter x > , called monomer 
activity, we define the following probability measure on the set S>g'- 



The normalising factor 



z g{x) := ^2 ■>' 



is called partition function of the model. Its natural logarithm log .2^ is 
called pressure. The expected value with respect to the measure hg,x is 
denoted by {-}q x , namely for any function / of the dimeric configuration 

(f)a, x --= E f{D)n G , x {D). 

Des> G 

Remark 1. The general monomer-dimer model on the graph G is obtained 
by assigning a monomeric weight x v > to each vertex v G V, a dimeric 
weight w e > to each edge e e E and considering the measure 

MG,x,w(-D) = Zg (^ >w ) lie eD IT* £^ G (D) ^ VD G QSq ■ 

In this paper we consider uniform monomeric and dimeric weights: x v = x, 
w e = w. Under this hypothesis one may assume without loss of generality 
w = 1, indeed it's easy to check that 



Z G (x,w) = W V\/*Z G (^,1). 



X 



Lemma 1. The pressure per particle admits the following bounds: 

1 \E\ 1 

logx < — \ogZ G {x) < \ogx+ — log(l + -3) . 
\V\ \V\ x z 



Proof. The lower bound is obtained considering only the empty dimeric 
configuration (i.e. a monomer on each vertex of the graph): 

Z G {x) > . 

The upper bound is obtained using the fact that any dimeric configuration 
made of d dimers is a (particular) set of d edges: 



\E\ \E\ 

Z G (x) = Card{D G @ G , \D\ = d} x^'™ < /~2 [ d ) 



x \V\-2d 



d=0 d=0 
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Remark 2. An important quantity of the model is the expectation of the 
fraction of vertices covered by monomers. A simple computation shows that 
it can be obtained from the pressure as: 

d \ogZ G (x) , \Jt G \ > 

We call this quantity monomer density. It is useful to introduce the following 
notation for the probability of having a monomer on a given vertex o G V: 

K x (G,o) :=(W G >G )a! e[0,l]. 

Now the monomer density can be rewritten as 

Two vertices u,v G V are neighbours in the graph G if there is an edge 
uv £ E connecting them: we write u ~ v . Denoting by E Q the set of edges 
which connect the vertex o£Vto one of its neighbours, we define the graph 
G-o:= (V ^o,E^E Q ). 

Following |9j we introduce a recursion relation for the probability TZ X (-) that 
will be extensively used in the sequel; this is a rewriting of the recursion 
relation for the partition function Z.{x) present in [2]. 

Lemma 2. The family of functions lZ x (G,o) fulfils the relation 

*- (G - o) = *» + S^£(g-°.«) < 5 » 

Proof. The dimeric configurations on G having a monomer on the vertex 
o coincide with the dimeric configurations on G — o. Instead the dimeric 
configurations on G having a dimer on the edge ov are in one-to-one corre- 
spondence with the dimeric configurations on G — o — v . Therefore 

1 V- „|.*c(D)| _ xZ G-o{x) 



;};'* 



Z G( X ) ot^ G Z G( X ) 

s.t. o^Mq(D) 



s.t. oe^#(j(-D) s.t.ooED 
= xZ G _ (x) + Y Z G-o-v(x) . 



\-* G {D)\ 



Hence one finds: 

v (r \ = xZ G _ (x) . Z G _ _ v {x) 

k x[ u, o, x + ^ Zg _ o _ v{x) \ L + ^v~o x Zq _ o{x) 



^ + x- 2 Y, v ~oK x (G-o,v)) 1 



■2 



x 2 + E v ~ KAG-o,v) ' 

□ 
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Iterating the recursion relation ([5]), one obtains immediately the squared 
recursion relation 

In the next lemma we allow the monomer activity to take complex values, 
precisely those of the open half-plane 

M + = {z G C|K(z) > 0}. 

This has no physical or probabilistic meaning, but it is a technique to obtain 
powerful results at real positive monomer activities by exploiting complex 
analysis. This lemma already appeared in [9j and in particular point ii can 
be seen also as a consequence of theorem 4.2 in [2]. 

Lemma 3. i. If z £ H+, then z" 1 K Z (G, o) G M + 

ii. The function z i— > lZ z (G,o) is analytic on M + 

Hi. If ze H+, then \TZ z (G,o)\ < \z\/St(z) 

Proof. Note that H + is closed with respect to the operations w h-> w~ l and 
(tOi,to 2 ) H> wi + w 2 . 

[i, ii] Proceed by induction on the number N = \ V\ of vertices of the graph 
G. For N = 1 the graph G coincides with its vertex o, hence 7Z Z (G, o) = 1. 
Therefore for 2 6 H+, aT 1 o) = z" 1 G M+ and K z (G,o) = 1 is 

obviously an analytic function of z. 

Suppose now the statements i and ii hold for any graph of N — 1 vertices 
and prove them for the graph G of iV vertices. By lemma [2j 



z 2 



By inductive hypothesis, for z G H + and for every v ~ o, z 1 TZ Z (G — o,v) G 
H_|_ and TZ Z (G — o,v) is an analytic function of z. Therefore 

z + Et>~o z ~ 1 T^z(G — o, v ) G H + (in particular it is 7^ 0) 

so that z~ 1 TZ z (G, o) G ELl and 1Z Z (G, o) is an analytic function of z (as it is 
the quotient of non-zero analytic functions). 

[iii] Use lemma [51 then apply the elementary inequality \z + w\ > 3i(z + w) 
and conclude using point i: 

\ n -AG, o)\ = \— r = -j- - - - I < 



2 + £„~ * _1 ^(c - M 1 - »( z ) + Ev-oH^ 1 - °. «)) 

S v ' 

\z\ >0 



In the graph G, given o G V and r G N, we denote by [G, o] r the ball of 
center o and radius r, that is the (connected) subgraph of G induced by the 
vertices at graph-distance < r from the origin o. 

A tree is a connected graph with no cycles. If the graph G is locally a tree 
near the vertex o, the next lemma allows to bound the operator lZ x {-,o) 
from above/below by cutting away the "non-tree" part of G at even/odd 
distance from o. 

Lemma 4 (Correlation inequalities on a locally tree-like graph). 

// [G, o]2r is a tree, then 1Z X {G, o) < 1Z X ([G, o]2 r , o) . 

If [G,o]2r+i is a tree, then TZ x {G,o) > TZ X ([G, o]2 r +i, o) . 

Proof. Proceed by induction on the distance r G N from the origin o. 
For r = 0, the graph [G, o]q is the isolate vertex o hence 

K x (G,o) <l = K x ([G,o\ ,o). 

Assume now the result holds for 2r and prove it for 2r + 1 (with r > 0). 
Suppose [G,o]2r+i is a tree. Note that [G, o]2 r +i — o = \_\ V ^ [G — o,v]2 r , 
where [G — o, v]2 r is a tree for every v ~ o. 

As in general 7Z x (H,v) depends only on the connected component of the 
graph H which contains the vertex v, it follows: 

Kx([G, o] 2r +i- o,v) = K X ([G - o,v] 2r ,v) . 

And by the induction hypothesis 

U X (G - o, v) < 1Z X ([G - o,v] 2r ,v) . 

Then using lemma [2] two times one obtains: 

x 2 x 2 



K x (G,o) = „ „ - - > 



x 2 + Ei~ K x (G-o,i) ~ x 2 + J2^ K x ([G-o,i]2r,i) 

TZ x ([G,o]2r+l,o) . 



x 2 



x2 + Ei~o^([ G '°]2r+l- 0,l) 

Induction from 2r — 1 to 2r (with r > 1) is done analogously. □ 



3 The model on a Galton- Watson tree 

Definition 3. As already said, a tree T is a connected graph with no cycles. 
A rooted tree is a tree T together with the choice of a vertex o, the root. 
This choice induces an order relation on the vertex set of T: the vertices 
which are neighbours of the root o form the 1 st generation, the vertices 
different from o and neighbours of a vertex in the 1 st generation compose 
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the 2 generation, and so on. Given a vertex v, its sons (or its offspring) 
are the vertices in the following generation which are neighbours of v. 
For r S N we denote T{r) the sub-tree of T induced by the vertices in the 
first r generations, namely T(r) = [T,o] r . The tree T is locally finite if the 
T(r)'s are finite graphs for every reN, 

The next proposition describes the behaviour of our model on any finite 
tree. It is an easy consequence of lemma [H 

Proposition 1. Let T be a locally finite tree rooted at o. Consider the 
monomer- dimer model on the finite sub-trees T(r), r G N. Then: 

i. r i — y lZ x (T(2r),o) is monotonically decreasing 

ii. r i — y lZ x (T(2r + 1), o) is monotonically increasing 
Hi. K x (T(2r), o) > ll x {T{2s + 1), o) Vr, s £ N 

Proof. Let r, s G N. 

[i] Consider the graph T(2r + 2). Cutting at distance 2r from o, one obtains 
[T(2r + 2), 6\2r = T{2r) which is a tree. Hence by lemmaH] 

K x {T(2r + 2),o) <K x (T(2r),o) . 

[ii] Consider the graph T(2r + 3). Cutting at distance 2r + 1 from o, one 
obtains [T(2r + 3),o]2r+i = T{2r + 1) which is a tree. Hence by lemma[3] 

K x (T{2r + 3), o) > 7^(T(2r + 1), o) . 

[iii] Consider the graph T(2r+1). Cutting at distance 2r from o, one obtains 
[T(2r + 1), o]2r = T(2r) which is a tree. Hence by lemmaS] 

K x (T(2r + l),o) <Tl x (T(2r),o). 

Now if r < s, combining point i. and this third inequality, one finds 

K x (T(2r), o) > K x (T(2s), o) > K x (T(2s + 1), o) ; 

while if s < r, combining point ii and the third inequality, one finds 

ft x (T(2s + l),o) <7^(T(2r + l),o) < TZ x (T{2r), o) . □ 

As a consequence of proposition Q] we obtain that on any locally finite 
rooted tree there exist lim r _>. 00 1Z x (T(2r), o) , limr^oo lZ x {T(2r + l),o) and 
moreover 

< lim K x (T(2r + l),o) = sup ^(T^r + 1), o) < 

< inf 7e x .(T(2r),o) = lim K x (T(2r),o) < 1. 

rgN r— >-oo 
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A natural question is if these two limits coincide or not. In the next 
proposition we prove that they are analytic functions of the monomer activ- 
ity x, so that it suffices to show that they coincide on a set of x's admitting 
a limit point to conclude that they coincide for all x > 0. We first state the 
following lemma of general usefulness. 

Lemma 5. Let (/ n )neN be a sequence of complex analytic functions onU C 
C open. Suppose that 

• for every compact K C U there exists a constant Ck < oo such that 

sn V \f n (z)\<C K VnGN; 

z&K 

• there exist Uq C U admitting a limit point and a function f on Uq such 
that 

fn(z) ►/(*) VzGC/q. 

n— >oo 

Then f can be extended on U in such a way that 

f n (z) ►/(*) VzGU; 

n— >oo 

further the convergence is uniform on compact sets and f is analytic on U . 

Proof. By hypothesis (/ n )neN is a family of complex analytic functions on 
U, which is uniformly bounded on every compact subset K C U. Therefore 
by Montel's theorem (e.g. see theorems 2.1 p. 308 and 1.1 p. 156 in [13J), 
each sub-sequence (/n m )meN admits a further sub-sequence (/n mp )peN that 
uniformly converges on every compact subset K C U to an analytic function 
/W, where a = (n mp ) peN . 

On the other hand by the second hypothesis one already knows that 

Vz G Uo 3 lim f n (z) . 

n— >oo 

Thus by uniqueness of the limit, all the /^ CT ^'s coincide on Uq . Hence, as 
Uo admits a limit point in U, by uniqueness of analytic continuation all the 
/(^'s coincide on the whole U. Denoting / their common value, this entails 
that 

VzeU 3 lim f n {z) = f{z). □ 
n— >oo 

Proposition 2. Let T be a locally finite tree rooted at o. Consider the 
monomer-dimer model on the sub-trees T(r), r G N. Then the maps 

x ^ lim K x (T(2r),o) , x ^ lim K x (T(2r +1), o) 
are analytic on M + . 
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Proof. Set f r (z) := TZ z (T(2r),o) and g r (z) := K K (T(2r + l),o). By lemma[3] 
(/r)rGN is a family of complex analytic functions on H+, and it is uniformly 
bounded on every compact subset K C H + : 



sup|/ r (z)| < sup-—- < oo VrGN. 
zeK zeK Ji{z) 

On the other hand by proposition Q] one already knows that 

Vx>0 3 lim f r (x) . 

r— )-oo 

The result for (/ r )reN then follows by lemma [5j The same reasoning holds 
for the sequence (fly)reN- D 

Now we define an important class of random trees. We will prove that for 
these trees the previous limits on even and odd depth almost surely coincide 
at every monomer activity. 

Definition 4. Let P = (Pk)k<=N > P = (Pfc)fceN be two probability distribu- 
tions over N . A Galton- Watson tree T(P, p) is a random tree rooted at o 
and defined constructively as follows. 

Let A be a random variable with distribution P, let (ifr,i)r>i,i>l be i.i.d. 
random variables with distribution p and independent of A . 

1) Connect the root o to A offspring, which form the I s * generation 

2) Connect each node (r, i) in the r th generation to K r ^ offspring; all 
together these nodes form the (r + l) th generation 

Repeat recursively the second step for all r > 1 and obtain T(P,p). 
We denote T(P, p, r) the sub-tree of T(P, p) induced by the first r genera- 
tions. Note that T(P, p) is locally finite. 

A special case of Galton- Watson tree is when p = P, which we simply denote 
T{p) ■= T(p, p) and TO, r) := T{p, p,r) . 

I f instead the offspring distributions satisfy P := Ylfc=a k Pk < oo and 

Pk = = VfceN, 

we call T(P, p) a unimodular Galton- Watson tree. 

In the following when we consider a Galton- Watson tree we suppose it is 
defined on the probability space (Q, T, P) and we denote E[ • ] the expectation 
with respect to the measure P. It is important to notice that when the 
monomer-dimer model is studied on a random graph G, then the measure 
Pg,x is a random measure and therefore the probability 1Z X (G, o) is a random 
variable. 
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Theorem 1. Let T(p) be a Galton-Watson tree such that ~p := EfceN k Pk < 
oo . Consider the monomer- dimer model on the finite sub-trees T(p,r), r G 
N. Then almost surely for every x > 

3 lim K x (T(p,r),o) =: X x . 

r— >oo 

The random function x \— > X x is almost surely analytic on M + . 

The distribution of the random variable X x is the only solution supported in 

[0, 1] of the following fixed point distributional equation: 

,,2 



V 

* 2 + Ef=i Xi 

where (Xi)i^ are i.i.d. copies of X, K has distribution p, and K 

are independent. 

Proof. To ease the notation we drop the symbol p as T := T(p) and T(r) := 
T(p,r). By proposition [1] there exist the two limits 

X+ := lim KJT(2r),o) , X~ := lim K x (T(2r + 1), o) , 

r—>oo r—too 

moreover < X~ < X£ < 1 and by proposition [2] the functions x \-t X£ 

and x I— > X~ are analytic on R+. 

The theorem is obtained by the following lemmata. 

Lemma 6. Given x > 0, X£ and X~ are both solutions of the following 
fixed point distributional equation: 

^^(i + E^^ + ESi^)- 1 )- 1 , (8) 

where (Xjj)jj e N are i.i.d. copies of X, {Hi)i^i are i.i.d. with distribution p, 
K has distribution p, (^ij)i.jeN; (-Hj)ieN an d K are mutually independent. 

We will write u <— v to denote "u son of v in the rooted tree (7~, o)". 
We will indicate 7L(r) the sub-tree of T induced by the vertex u and its de- 
scendants until the r th generation (starting counting from u). Using lemma 
[2] and precisely equation ([6]) one finds, with the notations just introduced, 

K x (T(2r + 2), o) = (1 + {x 2 + E u *-v K x (T(2r + 2)-o-v, u)) 1 ) 1 

= (i + £„<_„ (* 2 + E„^^(r u (2r) ! u))- 1 )- 1 

= (1 + E*i (x 2 + Ef4^(^(2r), ))- 1 )- 1 , 

where (7i,j(2r))ij 6 N are i.i.d. copies of T(2r), independent of (iTj)i 6 N an d 
if. 

Now since 7Z x (T(2r), o) a ' s ' > X+, it holds also 

J - — >CX> 

7^(r(2r),o)-^X+, 
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and moreover, thanks to the mutual independence of (TZ x (Tij(2r), o)) { . 
(iJj)j 6 pj , K , by standard probability arguments^ 



( [K x (Tij(2 r ), o)) , (fli)ieN , ( (X+ ) iJeN , (iJ^ew , # ) > 



where (X^jjgN are i.i.d. copies of X+, independent of (H,i)i € ^ and if. 
Then for any bounded continuous function <f> : [0, 1] — > M. 

E[<t>(X+)\ = hm E[<f>(TZ x (T(2r + 2),o))} 

r— >oo 

(i + Eli (* 2 + Efi^^^r)^))- 1 )- 1 ) 
-((i + E^^ + ESi^J)- 1 ) -1 ) 



lim E 

r— >-oo 



E 



Namely X+ is a solution of distributional equation (jSJ) . 
In an analogous way it can be proven that also X~ is a solution of distribu- 
tional equation ([8]). 

Lemma 7. Almost surely for all x > = X~ . 

By proposition [Tj X% > XjT. By lemma [6] and X~ are both so- 
lutions of equation (|SJ). Therefore, taking ((iTj)j 6 ^ , if) independent of 

((Xj)ijeN, (ly^jeN), one obtains: 



E[\X+-X~\} = |E[X+]-E[X-]| = 

= ^[(l+E^^+Efii^sr 1 ) -1 ] + 



-r 1 



E 



-E[(i + E£i(* 2 + Ef=iX 

E£i(^ 2 + Ef=i Xjj)' 1 - E£i(^ 2 + Ef=! x+r* . 
(i + E£i(^ 2 + Efi^i-)- 1 )!! + E£i(^ 2 + Efi-v,,) ') 



e[(e£x 



E,-=i(^?,- - x i, 



-) 



(^ 2 + Ef=i^)(x 2 + Ef=i^ 

• (i + Ef =1 (* 2 + Efi ^jrr 1 (i + Ef =1 (* 2 + Efii ^-rr 1 ] I 
< 4 E[Ef=i Ef=i 1*5 - ] = K n\x+ - x-\ } , 



where the last equality is true by independence. 

If x > \fj> , the contraction coefficient is p 2 /x 4 < 1. Therefore for all x > yfp 
E[|X+-X"|]=0, i.e. X+ = X~a.s. 



Equivalence between convergence in distribution and convergence of the characteristic 
functions (e.g. see theorems 26.3 p. 349 and 29.4 p. 383 in [T3]) can be used. 



12 



As Q is countable it follows that 

(X+ = X~ Vie]v^,oo[nQ) a.s. 

Now remind that by proposition [2] X~ are analytic functions of x > 0. 
Hence, as Q is dense in R, this entails that 

(X+ = X~ Vi > 0) a.s. 

by uniqueness of the analytic continuation. 

As a consequence f 3 lim T ._ > . 00 TZ x (T(r), 6) = X£ = X~ Vx > 0) a.s. We call 
X x this random analytic function of x. 

Lemma 8. Given x > 0, the random variable X x , satisfying the distribu- 
tional equation {5^, satisfies also the distributional equation (T7p. 

Using lemma [2] and precisely equation ([5]), one finds 
K x (T(r + l),o) 



v 



x 2 + Z v ^ K x (T(r + l)-o,v) x 2 + T, v <- K x (Tv(r),v) 

J2. 



"^ 2 + E£i^(^(r),o) ' 

where {%{r))i^ are i.i.d. copies of T(r), independent of K. 

Now since H x (T(r), o) °" s ' > X x (by definition, which is possible thanks to 

r— >oo 

lemma ED, it holds also 



K x (T(r),o) -^X x , 

r— yoo 



and moreover, thanks to the independence of {lZ x (Ti{r), o)) igN , K . 



{{K x (Ti(r),o)). &v K) -^-> ((Xi) im ,K), 

where (Xi)^ are i.i.d. copies of X x , independent of K. 
Then for any bounded continuous function (j) : [0, 1] — > M 

EW,)] = ME[^(T(r + l) )0 ))] = r^ Ki + ^f^mr)^ 

= m tv$k y )] • 

Namely Xj, is a solution of distributional equation ([7]). 

Lemma 9. For a given x > 0, the distributional equation ^ has a unique 
solution supported in [0, 1] . 
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Let Y be a random variable taking values in [0, 1] and such that 

Y v x2 



* 2 + E£i Yt ' 

where (Yi)i^ are i.i.d. copies of Y, independent of K. Observe that: 



x 2 



j'i- 
P>l 



x 2 + Ef=i ^ 



< 



Y K x (T(0),o) 

Therefore there exist (l^')ieN i-i-d- copies of Y and (7£ X (T(0), o)j)j g pj i.i.d. 
copies of T2.a;(T(0), o) such that 

^/ < K x (T(0),o)i ViGN. 



Let K' ^ p independent of (Y/)j 6 N ; (^x(7"(0), o)j)j e N . Applying the func- 
tion ^7 , which is monotonically decreasing in each argument, to 

* 2 +E*i(-) 

each term of the previous inequality one finds 

x 2 x 2 
< 



^ + Ef = ' 1 ^(r(o),o) ! ^ + E,=' 1 ^/ 
p>ii »n 

ft*(T(l),o) Y 

Therefore there exist (1Z X (T(1), o)i)j £ N i-i-d. copies of 7£ X (T(1), o) and 
(y/')j g N i.i.d. copies of Y such that 

ft:r(T(l), o)i < Y(' ViGN. 

Let If" ~ p independent of (72. X (T(1), o)i)ieN ? (X/O^eN • Applying the func- 

2 

tion 2 + ^k"( , which is monotonically decreasing in each argument, to 
each term of the previous inequality one finds 

x 2 x 2 



x 2 + ZZY>' " x 2 + ZZK x (T(l),o)> 
Oil P>ll 
Y K x (T(2),o) 

Proceeding with this reasoning one obtains that for any r 6 N there exist 
K x {T{r),o)~ = K x (T(r),o), Y^ = Y such that 

K x (T(2r + l),o)~ < y^+i) and Y^ < 1l x (T(2r),o)~ 

Pi I ftl 

4- as r — > oo 4- 
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Since by lemma [7J X x = X x = X x a.s., it follows^ that Y = X x . □ 

Corollary 1. Let T(P, p) be a Galton- Watson tree such that ~p := YlkeN k Pk < 
oo . Consider the monomer- dimer model on the sub-trees T(P, p,r), r £ N. 
Then almost surely for every x > 

3 lim K x (T(P,p,r),o) =:Y X . 

r— too 

The random function is a.s. analytic on R+. 

The distribution of the random variable Y x is 



* 2 + Ef=i Xi 

where A has distribution P and is independent of > (^Q)ieN are i-i-d. 

copies of X, the distribution of X is the only solution supported in [0, 1] of 
the following fixed point distributional equation: 



z 2 + £,=i^ ' 

where K has distribution p and is independent of (Xj)j e nj- 

Proof. We drop the symbols P, p as T* := T(P,p) and T*{r) := T(P,p,r). 
Observe that T* — o = \_\ v4 _ T^ and the random trees (T^) v <r-o are i.i.d. 
Galton- Watson trees of the type T(p). Using lemma [2] 



Tl x {T*{r + l),o) 



x 2 x 2 



x 2 + Ev<-oK*<T*(r + 1) - o,v) x 2 + Ev<-oK*(TS(r),v) 



By theorem [T] for any v son of o, lim^oo 1Z x (Ty(r), o) almost surely ex- 
ists, it is analytic, and its distribution satisfies equation ([7]). Therefore 
lim r _ ) . 00 TZ x (T*(r), o) almost surely exists and is analytic, in fact 



lim TZ x {T*{r),o) 



x 2 



x2 + E«<- lim r ^oo Tlx (T*(r),v) 

v 



„2 



^ 2 + Ef=i^' 



where pQ)i e nj are i.i.d. copies of the solution supported in [0, 1] of equation 
([7]) , A has distribution P and is independent of (Xj)j G N • D 



2 A squeeze theorem for convergence in distribution holds: if X n < Y„, < X' n , 
Y n = Y^ = Y for all n £ N and X n X, X' n — X then Y = X. 

To prove it work with the CDFs: F x ,Jx) < Fy^(x) = Fy n {x) < F Xn (x) Vx g R, 
Fx' (x) > Fx (x) and Fx n (x) > Fx (x) for every x continuity point of Fx ■ Since 

n n— >-oo n — ^oo 

Fy = Fy n = Fy, by the classical squeeze theorem it follows that Fy(x) = F x (x) for every 
x continuity point of Fx. Now since Fx and Fy are right-continuous and the continuity 
points of Fx are dense in R, one concludes that Fy = Fx- 
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Corollary 2. In the hypothesis of corollary [Jl almost surely for every z 6 
3 lim K z {T{P,p,r),o) =:Y Z . 

1 — >oo 

The random function z \-± Y z is almost surely analytic on H + and the con- 
vergence is uniform on compact subsets ofW + . 

Proof. Set fr(z) := TZ Z (T(P, p,r),o). By lemma (/r)reN is a sequence of 
complex analytic functions on H+, uniformly bounded on compact subsets. 
On the other hand by corollary [T] (/ r )reN a - s - converges pointwise on K + . 
Then the result follows from lemma [5j □ 



4 The model on random graphs locally convergent 
to a Galton- Watson tree 

Let G n = (V n ,E n ), n e N be a sequence of random finite simple graphs, 
defined on the probability space (Q, F, P). 

We introduce now the main class of graphs studied in this paper. The idea 
is to fix a radius r and draw a vertex v uniformly at random from the 
graph G n : for n large enough wg want the ball [G^,, v] r to be a (truncated) 
Galton- Watson tree with arbitrary high probability. 

Definition 5. The random graphs sequence (G n ) n£ N locally converges to 
the unimodular Galton- Watson tree T(P, p) if for any r £ N and for any 
(T, o) finite rooted tree with at most r generations 

£l(([G B> «] r ,«)S(T )0 )) ^ P((T(P,p,r), ) -(!>)) 

1 n| i>ev„ 

Here = denotes the isomorphism relation between rooted graphs. 
Remark 3. The following statements are equivalent: 

i- (G n )neN locally converges to T(P,p) 

ii. a.s. for all r 6 N and (T, 6) finite rooted tree with at most r generations 

J2H([GnMr:V)^(T,0)) — -> P((T(P, /?, r), o) = (T, o)) 

iii. a.s. for all r 6 N and F bounded function, invariant under rooted 
graph isomorphisms, 

-^i V F([G n ,v] r ,v) l([G n ,u] r is atree) ► E[F(T(P, p, r), o)] 
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iv. a.s. for all r G N and (B,o) finite rooted graph with radius < r 

Vl(([G n ,u] r ,t;)^(B,o)) ► ¥((T(P,p,r),o)^(B,o)) 

vev n 

v. a.s. for all r G N and F bounded function, invariant under rooted 
graph isomorphisms, 

-i- V F([G n ,v]r,v) ► E[F(T(P,p,r),o)] 

Proof. Let (r), <^(r) be respectively the set of finite rooted graphs, trees 
with radius < r, considered up to isomorphism. It is important to note that 
they are countable sets. In particular let %(r), 5^(r) be respectively the 
set of finite rooted graphs, trees with radius < r and maximum degree < d, 
and observe that 

• ^d( r ) an d ^d( r ) are finite, indeed they contains only graphs with at 
most 1 + al + (d - l) 2 H h {d - l) r vertices 

• %(r) C % +1 (r) and 5^(r) C ,% +l {r) , 

• S?(r) = IUen^W and <T(r) = U deN W • 

We are interested in the two following probability measures on W(r) 

v r ,n(B,o) := £l(([G n ,i;] r ,t;)?*(B, )) V (B, o) € Sf (r) , 

i/ r (B, o) := P((T(P, p, r), o) (B, o)) V (B, o) G Sf (r) . 

Note that z^ r ,n is a random measure since it is an empirical average over the 
balls of the random graph G n . Fixed an elementary event uj G Q, we write 
v^ n for the corresponding deterministic measure. Note instead that v r is a 
deterministic measure supported on 3T(r). 

[i => ii] By hypothesis i for all r G N and (T, o) G ^(r) there exists a 
measurable set N r tT,o) suc h that P(AT r ; (t,o)) = an d 

< n (T, o) ► i/ r (T, o) Vw \ iV ri(T)0) . 

As UreN ^( r ) is countable, setting N := |J reN U(i»e^(r-) ^.(T.o) we obtain 
that P(iV) = and 

i? n (T,o) >v r {T,o) V(T,o)G^(r) Vr G N VwG^xiV. 

[ii =>• iii] By hypothesis there exists iV with F(N) = such that for all 
u G fi \ iV, r G N, (T, o) G ^(r) 

< n (T,o)— >i/ r (T,o). 
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Now let u G SI \ N, r G N and F : Sf(r) — > R bounded. Summing over 7d(r) 
which is finite, clearly: 



F(T,o)v» n (T,o) — -> ^ F(T, )i/ r (T, ). 

(T,o)e5i(r) (T,o)e5i(r) 

On the other hand the sum over the countable set «^'(r) := ^(r) \ J^(r) 
is: 

| F ^ °Xn{T,o) | < SUp |F| ^ <nC?»< 

(T,o)e^'(r) (T,o)6^'(r) 

sup|F|(l- £ < n (T,o)) — sup|F|(l- £ „ r (T >0 )) — > 0. 

(T,o)e5i(r) (T,o)e^(r) 

where the limit in n is done by hypothesis u and finiteness of ^(r), while 
the limit in <i is done by monotone convergence. Similarly one finds that: 

| Yl F(T,o)v r (T,o)\ < sup|F| MT,o) < 

(T,o)e^(r) (T,o)e^(r) 

sup|F|(l- V i/ r (T,o)) ► 0. 

' - , cf-s-oo 
(T,o)e£T d (r) 

These tree facts prove (using triangular inequality and lim sup) that 

Y F(T,o)^ n (T,o) Yl F(T,o)u r (T,o). 

(T,o)ef(r) n °° (T,o)e&(r) 

[iii =>- iv] By hypothesis in there exists N with P(iV) = such that for all 
u; G ft \ iV, r G N, F : &(r) -»• R bounded 

£ F(T, )< n (T, ) — > ^ F(7>K(T, ). 

(T,o)6^(r) " °° (T,o)€^(r) 

Let w G \ TV, r G N, (T, o) G ^(r). Taking F(-) = 1(- 9* (T, o)) , clearly 

< n (T, ) — - -> i/ r (T, ). 

Let instead (-B,o) G Sf(r) \ ^(r). Clearly u r (B,o) = and on the other 
hand, taking F = 1, 

< n (B,o) < 1- E <„(r,o) — ->• 1- ^ i/ r (T,o) = 0. 

(r,o)6^(r) (T,o)e^(r) 

[iv =>■ v] This proof is very similar to ii => m. By hypothesis if, there exists 
iV with P(AT) = such that for all u G fi \ TV, r G N, (5, o) G Sf(r) 



< n (S,o) — -> i/ r (5,o) 
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Now let oj G f2 \ N, r G N and F : t S{r') — > R bounded. Summing over %(r) 
which is finite, clearly: 

]T F(B j0 )z£ n (B,o) ► J] F(B,o)u r (B,o). 

(B,a)e%(r) (B,o)e%(r) 
On the other hand the sum over the countable set &'Ar) '■= @(r) \ is: 

| F(B,o)v? >n (B,o)\<sup\F\ <n( B >°) = 

sup \F\(1- V < n (S,o)) ► sup|F|(l- V v r (B,o)) > 0, 

(B,o)e%(r) (B,o)£%(r) 

where the limit in n is done by hypothesis iv and finiteness of %(r), while 
the limit in d is done by monotone convergence. Similarly one finds that: 

| F(B,o)u r (B,o) | < sup \F\ ^2 M B ,o) < 

{B,o)&%{r) (B,o)&fi(r) 

sup|F|(l- V uJB,o)) ► 0. 

(B,o)eSf d (r) 

These tree facts prove (using triangular inequality and lim sup) that 
V F(B,o)v? >n (B,o) ► V F(B,o)u r (B,o). 

(B,a)&f(r) (B,a)&f(r) 

[v i] By hypothesis u there exists N with P(iV) = such that for all 
uGO\JV,rGN,F: (r) -> R bounded 

^ F(B j0 )i£ n (B,o) ► Yl F(B,o)u r (B,o). 

(B,a)&f(r) (B,a)&f(r) 

Let co G tt n iV, r G N, (T,o) G ^(r). Taking F(-) = 1(- = (r,o)) , clearly 

< n (T,o) ► ^(T,o). □ 

Observe that local convergence of random graphs (G n ) ne ^ to the ran- 
dom tree T(P,p) is, in measure theory language, a. s.— weak convergence of 
random measures (^V,n)neN to the measure v r for all r G N. Prom this point 
of view remark [3] gives different characterisations of the weak convergence of 
measures, valid in general for measures defined on a discrete countable set (in 
particular the equivalences ii 44> Hi and iv 44> v can be seen as consequences 
of the Portmanteau theorem, e.g. see theorem 2.1 p. 16 in [15J). 
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Remark 4. In a graph G the degree of a vertex v, denoted deg G (i>), is the 
number of neighbours of v. If (G n ) ne ^ locally converges to T(P,p), then P 
is the empirical degree distribution of G n in the limit n — > oo. 
Indeed the degree is a local function (deg G (v) = degiQ ^^v)) and clearly an 
indicator function is bounded, hence by remark [3] a.s. for every k G N 

2 l(deg G » = k) —-^ P(deg r(Pp) ( ) = k) = P k . 

Definition 6. The random graphs sequence (G n ) ne N is uniformly sparse if 
lim limsup — — V deg Gn (w) l(deg Gn (u) > I) = a.s. 

Remark 5. If (G n ) ng N is uniformly sparse and locally convergent to T(P, p), 
then 

\EJ 1 ^ 



-)■ - P a.s. 



To prove it we write two times the number of edges as the sum of all vertices' 
degrees 

2 w\ = W~\ Y, de ^M- 

Then we fix I G N and we split the right-hand sum in two parts, concerning 
respectively smaller and grater than I degrees . To the first part we can 
apply the local convergence hypothesis (remark [3]): 

-L Y, de gG» l(deg G » < E[deg r(Pp) (o) l(deg r(Pp) (o) < I)] 

\v n \ veVn 

>E[deg r(Pp) (o)] =P. 

To the second part we apply the uniform sparsity hypothesis: 

lim limsup — Y] deg Gn (v) l(deg Gn (v) > 1+1) = a.s. 

l—>co n ->oo \V n \ *— ~ 
v£V n 

Example 1. An Erdds-Renyi random graph G n is a graph with n vertices, 
where each pair of vertices is linked by an edge independently with proba- 
bility c/n. The sequence (G n ) ne ^ is uniformly sparse and locally converges 
to the unimodular Galton- Watson tree T(P,p) with P = p = Poisson(c). 
For proof and further examples see |16t [7j . 

The next theorem is the main result of this paper. 
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Theorem 2. Let (G n ) ng pj be a sequence of finite random graphs, which: 

i. is locally convergent to the unimodular Galton- Watson tree T(P,p); 

ii. has asymptotic degree distribution P with finite second moment ( equiv- 
alently ~p < oo). 

Consider the monomer- dimer model on the graphs G n , n 6 N. Then almost 
surely for all x > the monomer density 

vev n 

The function x t— >■ E[Y X ] is analytic on M + . 

The random variable Y x is defined in corollary {J\ that is its distribution is: 

Y V * 



* 2 + £f=i^ 



where A has distribution P and is independent of {Xj)i^ , (Xj)jgN are i.i.d. 
copies of X, the distribution of X is the only solution supported in [0, 1] of 
the following fixed point distributional equation: 



* 2 + E£=i Xt ' 
where K has distribution p and is independent of (Xi)i^. 
Proof Set T* := T(P,p) and T*(r) := T(P,p,r). 

Let r G N and v G V n . If [G n ,v]2 r +i is a tree, then lemma H] permits to 
localize the problem: 



U x (G n ,v) t([G n ,v] 2r +i is a tree) 



< K x ([G n ,v] 2r ,v) l([G n ,v] 2r +i isatree) 
> TZ x ([G n ,v} 2r +i,v) l([G n ,v] 2r+ i isatree) 



Now work with the right-hand bounds and take the averages over a uniformly 
chosen vertex v. First let n — > oo using the hypothesis of local convergence 
(see remark [3|) and then let r — > oo using the results on Galton- Watson trees 
(corollary [1]) and dominated convergence: almost surely for all x > 

— r TZ x ([G n ,v] 2r ,v) t([G n ,v] 2r+ i is a tree) > 

Vn * ' ' 71— 5-CXD 



E[K x (T*(2r),o)] \ E[Y X ] 
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and similarly 

Tw~\ K x ([G n ,v]2r+i,v) l([G n ,v} 2r +i is a tree) 
lVnl vev n 

E[K x (T*(2r + l),o)] / E[Y X ]. 

r— >oo 

On the other hand observe that a.s. for all x > 



7777 K x (G n ,v) - -— ^2 K x (G n ,v) t([G n ,v} 2r +i is a tree) I < 
-Vfl- l([G n ,ul 2 r+l is a tree)) ► 1 -P(T*(2r+ 1) is a tree) = 0. 



Therefore one finds that almost surely for all x > 

limsup— V TZ x {G n ,v) < E[Y X ] ; liminf — , V K x {G n ,v) > E[Y X ) . 

Namely there exists 

lim — 1 V TZ x (G n ,v) = E[Y X ) Vx>0 a.s. 

Remembering remark [2] and in particular the identity (JH) the proof is con- 
cluded, except for the analyticity of x 1— > E{Y X ] which will follow from the 
next corollary. □ 



Corollary 3. In the hypothesis of theorem^ almost surely for all z € 

v€V„ 

where the random variable Y z is defined in corollary^ 

The function z 1— > E[Y Z ] is analytic on H + and the convergence is uniform 
on compact subsets of M + . 

As a consequence almost surely for all k > 1 and z E H_|_ 
d fc log Z Gn {z) > d k E[Y Z ) 

dz k \V n \ n->oo dz k z 

Proof. By lemma [3] (£G„)neN is a sequence of complex analytic functions on 
HL|_, which is uniformly bounded on compact subsets K C EI+: 

sup|eG„(^)| < ttt-t Y] sup \TZ z (G n ,o)\ < sup^-- < 00 VnGN. 
z€K \V n \ z£K zeK Ji(z) 
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On the other hand by theorem [2] (s Gn (ic)) n eN a.s. converges pointwise on 
K_j- to E[YJ. Then lemma [5] applies: E[Y^] is analytic in z 6 H + and a.s. 

£g n (z) t- E[Y 2 ] uniformly in z € K for every compact K C H_|_. 

This entails also convergence of derivatives (e.g. see theorem 1.2 p. 157 in 
EL31). □ 



The existence and analyticity of the monomer density in the thermody- 
namic limit entails the same properties for the pressure per particle. Only 
the additional assumption of uniform sparsity is required. 

Corollary 4. Let {G n ) n ^ be a sequence of random graphs, which: 

i. is locally convergent to the unimodular Galton- Watson tree T(P,p); 

ii. has asymptotic degree distribution P with finite second moment; 
Hi. is uniformly sparse. 

Then almost surely for every x > 

e[y] 



\v, 



i r 

— logZ Gn (x) p{a) + / 



dt 



where a > is arbitrary, p{a) = linin^oo ^r-r log Z Gn (a) a.s., and Y t is the 

random variable defined in theorem^ 

The function x i— > p(a) + J x dt is analytic on R + . 

Proof. From theorem^ using the fundamental theorem of calculus and dom- 
inated convergence, it follows immediately that a.s. for every x > 0, a > 

logZ Gn (x) log Z Gn (a) _ [* d logZ Gn (t)_ dt > rMM d i ( 9 ) 



\v n \ \v n \ J a dt \v n 

By theorem[2]the function x \— > E[Y X ] is analytic on IR_|_, therefore the integral 

function x i— >• J x dt is analytic on M + too. 

To conclude it remains to prove that almost surely for all x > 

3 lim lQg ^ (x) . 

n— too \Vn\ 



Use the bounds for the pressure of lemma [T] to estimate 
logZ Gn (x) log Z Gn (a) 



\V„\ \V„\ I > l2i^M _l 0ga _^l 0g( i + J, 



(10) 
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Put together ©, (ITUL remind |P n |/|y n | — — > P/2 and obtain that a.s. for 

ra— )-oo 

all x > 



iiminf log ^w >io go+ rsa dt , 

n->oo \V n \ J a t 

log Z G (x) , P , . 1 x r EM , 

limsup , TP ; < loga+- log(l + — ) + / -J^dt. 

Therefore a.s. for all x > 

< hmsup l0g ^7 (x) -liminf l0g ^^ < £ log(l + 1) .0, 

n->oo \V n \ n->oo |l/ n | 2 GT a-^oo 

which entails existence of lim^^oo Iog f*?*? ^ and completes the proof. □ 

Corollary 5. In the hypothesis of corollary [^} if P > 0, almost surely the 
pressure 

log Z Gn 



lim 



n— >oo K, 



is an analytic function of the monomer density 

lim e Gn . 

n— >oo 

Proof. Set p n := °| V J" , p := lim n ^ 00 p n and e n := e Gn , e := ]jm n -+ O0 £ n . 
By theorem [2] and corollary U on an event of probability 1 the monomer 
density e and the pressure p are analytic functions of the monomer activity 
x > 0. Now a direct computation shows that 

8e n . < \^G n \ 2 >G n ,x ~ < \Jfc n \ >G n ,x . n 

But a more precise lower bound is provided by theorems 7.3 and 7.6 in [2J: 

x^(x) > ^x 2 (l-e n {x)) 2 and 1 - e n (x) > -^-^ , 
hence 

de n Ax 2 \E n \ 2x 2 - 

(x 2 + 2) 2 |K| (z 2 + 2) 2 

By corollary [3] it follows: 

de. N 2x 2 - 

ax (x z + 2)^ 

Thus e is an analytic function of x with non-zero derivative, so that it is 
invertible and its inverse is analytic (e.g. see theorem 6.1 p. 76 of |13|). In 
other words x can be seen as an analytic function of e. Since the composition 
of analytic functions is analytic, it is proved that p is an analytic function 
off-. □ 
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5 Upper and lower bounds 



To conclude we consider the particular case when the graphs sequence 
(G n ) ng N locally converges to T(P, p) with P = p = Poisson(2) (e.g. this is 
the case of G n Erdos-Renyi with c = 2), and we show an approximate plot 
of the monomer density e(x) := lim eg„(x)- 




Fig. 2: Upper and lower bounds for the monomer density e on the Erdos-Renyi 
with c = 2 versus the monomer activity x (squares, triangles, diamonds, circles). 
The monomer density on the binary tree (continuous line) and the complete graph 
(dashed line) versus the monomer activity, see 



We describe briefly how to obtain it. The distributional recursion X — 
x 2 / (x 2 + YliLi Xi) with K ~ P = Poisson(2) is iterated a finite number r of 
times with initial values Xi = 1. The obtained random variable X(r) rep- 
resents the monomer density on a truncated Galton- Watson tree T(P, P, r) 
(lemma[2]). If X is the fixed point of the equation, we know that X(2r) \ X, 
X(2r + 1) /* X as r — > oo (proposition [H theorem [1]) and that K[X] is the 
asymptotic monomer density on (G n ) ne js} (theorem [2]) . 
For values of x = 0.01, 0.1, 0.2, . . . , 2, the random variables X(r), r = 
3, 4, 5, 6 are simulated numerically 10000 times and an empirical mean is 
done in order to approximate E[X(r)]. The results are plotted as circles, 
squares, diamonds, triangles connected by straight lines. 
The dot at 0.216074 on the vertical axes corresponds to the exact value of the 
monomer density when the monomer activity x — > 0, supplied by the Karp- 
Sipser formula [10] or by its extension due to Bordenave, Lelarge, Salez [9j. 
Therefore the graph of the monomer density x i— > E[X] = lim eg„ (x) starts 

n— ¥oo 

from (0, 0.216074) and lays between the diamonds and triangles curves. 
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